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' | ' Abstract 

This paper is divided into two parts. In the first part we analyze the consequences, 
for the LHC, of gauge and third family Yukawa coupling unification with a particular 
set of boundary conditions defined at the GUT scale. We perform a global x 2 anal- 
ysis including the observables M\y, Mz,Gf,o;~^, a s (Mz), M t ,mb(mb), M T , BR(B — > 
X a i), BR(B S ->n+n-) and M h . The fit is performed in the MSSM in terms of 9 GUT 
CN scale parameters, while tan (3 and /x are fixed at the weak scale. Good fits suggest an up- 

per bound on the gluino mass, M g < 2 TeV. This constraint comes predominantly from 
q fitting the bottom quark and Higgs masses (assuming a 125 GeV Higgs). Gluinos should 

^ be visible at the LHC in the 14 TeV run but they cannot be described by the typical 

t-h simplified models. This is because the branching ratios for g — >■ tt Xi,2-> ^ Xi,2> ^ X121 

bi X121 9 X1234 are comparable. Stops and sbottoms may also be visible. Charginos 
and neutralinos can be light with the LSP predominantly binodike. In the second part 
.£h of the paper we analyze a complete three family model and discuss the quality of the 

global x 2 fits and the differences between the third family analysis and the full three 
family analysis for overlapping observables. We note that the light Higgs in our model 
couples to matter like the Standard Model Higgs. Any deviation from this would rule 
out this model. 



1 Introduction 



Gauge coupling unification in supersymmetric grand unified theories (SUSY GUTs) [1-6] provides an 
experimental hint for low energy SUSY. However, it does not significantly constrain the spectrum of 
supersymmetric particles. On the other hand, it has been observed that Yukawa coupling unification 
for the third generation of quarks and leptons in models, such as SO(10) or SU(4) c xSU(2) l xSU(2)#, 
can place significant constraints on the SUSY spectrum in order to fit the top, bottom and tau 
masses [7-11]. These constraints depend on the particular boundary conditions for sparticle masses 
chosen at the GUT scale (see for example, [9, 12, 13], which consider different GUT scale boundary 
conditions). In light of the present success of the LHC with the possible observation of the Higgs 
boson with mass of order 125 GeV and significant lower bounds on gluino and squark masses, it 
is a perfect time to review the viability of the constraints on the sparticle spectrum resulting from 
gauge and third generation Yukawa coupling unification. 1 This is what we do in this paper. In 
part one of the paper, we perform a global x 2 analysis assuming SO(10) boundary conditions for 
sparticle masses and non-universal Higgs masses, which we have called "just so Higgs splitting." We 
fit the observables, M w , M z , G F , a~^, a s (M z ), M t ,m b (m b ), M T , BR(B -> X s >y),BR(B s -)■ 
and Mh in terms of 11 arbitrary parameters. These fits then place significant constraints on the 
gluino mass. 

In the second part of the paper we study a complete three family model of quark and lepton Yukawa 
couplings at the GUT scale [16, 17] which is based on an SO(10) GUT with &D 3 x [U(l) xZ 2 x Z 3 ] 
family symmetry. This model was shown to give good fits to precision electroweak data, including 
quark, charged lepton and neutrino masses and mixing angles (see most recently the global x 2 
analysis in [18]). In light of the observation of sin 2 #13 it is again a perfect time to re-analyze 
this model. We are also able to compare the third family Yukawa unification analysis with the 
three family analysis which now includes hierarchical Yukawa matrices with unification of the (3,3) 
element of the Yukawa matrices. Hence, off-diagonal elements in the Yukawa matrices give small 
corrections to exact Yukawa unification. 

The paper is organized as follows. In Section 2, we present the SO(10) model. In Section 3, we 
present the procedure used in the paper for analyzing the model. In Section 4, we consider a model 
with gauge coupling unification and only the Yukawa couplings for the third family, which are 
assumed to unify at the GUT scale. We perform a global x 2 analysis fitting the relevant low energy 
observables. In Section 5, we extend the analysis to all three families of quarks and leptons using a 
particular SO(10) GUT model. In this case, we look for the minimum values of x 2 f° r h ve different 
choices of the universal squark and slepton mass, mig, defined at the GUT scale, Mq. Finally, the 
summary and conclusions are given in Section 6. 

For other analyses in this direction, see [14, 15]. 
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2 The Model 



2.1 Third family model 

Fermion masses and quark mixing angles are manifestly hierarchical. The simplest way to describe 
this hierarchy is with Yukawa matrices which are also hierarchical. Moreover the most natural way 
to obtain the hierarchy is in terms of effective higher dimension operators of the form 

45 

W D A 16 3 10 16 3 + 16 3 10 — 16 2 + • • • • (1) 

M 

This version of SO(10) models has the nice features that it only requires small representations of 
SO(10), has many predictions and can, in principle, find an UV completion in string theory. The 
only renormalizable term in W is A I63 10 I63 which gives Yukawa coupling unification 

A = At = \b = A r = X Ut (2) 

at Mqut- Note, one cannot predict the top mass due to large SUSY threshold corrections to the 
bottom and tau masses, as shown in [19-21]. These corrections are of the form 

azixM-g tan/3 A? fi A t tan/3 

drub/mi, oc % 1 k h log corrections. (6) 

m~ mf 
b t 

So instead we use Yukawa unification to predict the soft SUSY breaking masses. In order to fit the 
data, we need 

5m b /m b ~ -2%. (4) 

We take /i, M g > 0, thus we need fi, A t < 0. For a short list of references on this subject, see 
[7-11, 14, 22-24]. 

Given the following GUT scale boundary conditions, namely universal squark and slepton masses, 
mi6, universal cubic scalar parameter, Aq, universal gaugino masses, M 1 / 2 , and non-universal Higgs 
masses [NUHM] or "just so" Higgs splitting, m,H u , m H d or m 2 H = m 2 [l — (+)A^ H ], we find that 
fitting the top, bottom and tau mass forces us into the region of SUSY breaking parameter space 
with 

A m — 2mi6, mio ~ y/2 mie, mie > few TeV, fi,My 2 <C mi 6 ; (5) 

and, finally, 

tan/3«50. (6) 

In addition, radiative electroweak symmetry breaking requires ~ 13%, with roughly half of 

this coming naturally from the renormalization group running of neutrino Yukawa couplings from 
Mq to M Nt ~ 10 13 GeV [9]. 

It is very interesting that the above region in SUSY parameter space results in an inverted scalar 
mass hierarchy at the weak scale with the third family scalars significantly lighter than the first 
two families [25]. This has the nice property of suppressing flavor changing neutral current and CP 
violating processes. These results depend solely on SO(10) Yukawa unification for the third family. 
In order to demonstrate this, we perform a separate analysis with only third family observables 
(Section 4) and then a complete three family analysis (Section 5). 2 



2 The large Yukawa coupling for the third family is the driving force for the inverted scalar mass hierarchy. However, 
the particular boundary conditions of Eq. (5) were shown to maximize the effect. 
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2.2 Full Three Family Model 



We now consider a complete three family SO(10) model for fermion masses and mixing, including 
neutrinos [16-18]. The model also includes a D% x [U(l) x Z2 x Z3] family symmetry which is 
necessary to obtain a predictive theory of fermion masses by reducing the number of arbitrary 
parameters in the Yukawa matrices. Consider the superpotential generating the effective fermion 
Yukawa couplings: 



W ch . 

fermions — 

A I63 10 I63 + 16 a 10 Xa + Xa {M % X a + 45 



™ 16 3 + 45 ^ 16 Q + A 16 a ) (7) 
M M 



where 45 is an SO(10) adjoint field which is assumed to obtain a VEV in the B - L direction, 
M x is a linear combination of an SO(10) singlet and adjoint, and the index a = 1,2. Its VEV 
M x = Mq(1 + aX + f3Y) gives mass to Froggatt-Nielsen states [26]. Here, X and Y are elements of 
the Lie algebra of SO(10) with X in the direction of the U(l) which commutes with SU(5) and Y 
the standard weak hypercharge, and a , j3 are arbitrary constants which are fit to the data. M is 
an SO(10) invariant mass scale which in principle could be obtained by integrating out additional 
Froggatt-Nielsen states. Note that both Mq and M are assumed to be above the GUT scale. (fi a , 
4> a , A are SO (10) singlet "flavon" fields, A is a non-trivial one dimensional representation under 
Z?3, and Xa, Xa are a pair of Froggatt-Nielsen states transforming as a 16 and 16 under SO(10). 
The so-called flavon fields are assumed to obtain VEVs of the form 







(8) 



After integrating out the Froggatt-Nielsen states one obtains the effective fermion mass operators 
in Fig. 1. 
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16, 



16, 




-X- 



(a) Renormalizable mass term for third family that (b) Effective operator that generates the (2, 2) element 
gives rise to the (3, 3) element of the Yukawa matrix, of the Yukawa matrix. 




M, 



16 b 



16c 



(c) Effective operators that generate the off-diagonal (d) Effective operators that generate the off-diagonal 
Yukawa couplings: (b, c) = (3,2), (2,3), (3,1) or (1,3) Yukawa couplings: (b,c) — (2,1) or (1,2). 

Figure 1: The effective fermion mass operators obtained after integrating out the Froggatt-Nielsen massive states. 
Here, a runs from 1 to 2. 



Inserting the flavon VEVs, one then obtains Yukawa matrices for up-quarks, down-quarks, charged 
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leptons and neutrinos given by 



Y u 




Y d = -e> e -ea A (9) 



with 



-3 e £ cr 



£ = fa/ fa, i a 02/M, (10) 
e oc fa/M, e' ~ (A/M ), 



1 + a 
1 -3a' 



and 




F„ = e' u 3 e w few A (11) 



with co = 2a /{2a — 1) and a Dirac neutrino mass matrix given by 

m v = Y u ^= sin /3. (12) 
\J2 

From Eq. (9) and Eq. (11) one can see that the flavor hierarchies in the Yukawa couplings are 
encoded in terms of the four complex parameters p, a, e, £ and the additional real ones e, e', A. These 
matrices contain 7 real parameters and 4 arbitrary phases. Note, the superpotential (Eq. (7)) has 
many arbitrary parameters. However, at the end of the day the effective Yukawa matrices have 
much fewer parameters. This is good, because we then obtain a very predictive theory. Also, the 
quark mass matrices accommodate the Georgi-Jarlskog mechanism, such that m^/m e « 9m s /md. 
This is a result of the 45 VEV in the B — L direction. 

We then add 3 real Majorana mass parameters for the neutrino see-saw mechanism. The anti- 
neutrinos obtain GUT scale masses by mixing with three SO(10) singlets (N a for a = 1,2 and 
A^) transforming as a D% doublet and singlet, respectively. The full superpotential is given by 
W = W ch . 

fcrmions + ^neutrino witn 

Neutrino = 16 (A 2 N a 16 a + A 3 N 3 I63) + \ (S a N a N a + S 3 N 3 A3) (13) 

where the fields S a , S3 are additional flavon fields whose VEVs provide Majorana masses for the 
states A a , A3. We assume 16 obtains a VEV, viq, in the right-handed neutrino direction, and 
(S a ) = M a for a = 1, 2 and (S3) = M3. The effective neutrino mass terms are given by 

W = vm v u + uVN+]^NM N N (14) 
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Sector 


Third Family Analysis 


# 


Full three family Analysis 


# 


gauge 


a G , M G , e 3 


3 


a G , M G , e 3 


3 


SUSY (GUT scale) 


""116, My 2 , A , rn Hu , rn Hd 


5 


mi 6 , M 1/2 , M, m H u , m Hd 


5 


textures 


A 


1 


e, e', A, p, a, e, £ 


11 


neutrino 







M Rl , M R2 , Mr 3 


3 


SUSY (EW scale) 


tan /3, p 


2 


tan/5, p 


2 


Total # 




11 




24 



Table 1: The model is defined by three gauge parameters, a G ,M G ,e3; one large Yukawa coupling, A; 5 SUSY 
parameters defined at the GUT scale, mie (universal scalar mass for squarks and sleptons), Mi/ 2 (universal gaugino 
mass), m,H u , m H d (up and down Higgs masses), and Aq (universal trilinear scalar coupling); p, tan/? obtained at 
the weak scale by consistent electroweak symmetry breaking. The full three family model has additional off-diagonal 
Yukawa couplings, and includes 3 right-handed neutrino masses. 



with 

/ A 2 

V = v 16 A 2 | , M N = diag(Mi, M 2 , M 3 ) (15) 
V A 3 

all assumed to be real. Finally, upon integrating out the heavy Majorana neutrinos we obtain the 
3x3 Majorana mass matrix for the light neutrinos in the lepton flavor basis given by 

M = Uj m u M R l ml U e , (16) 

where the effective right-handed neutrino Majorana mass matrix is given by: 

M R = V M N l V T = dmg(M Rl ,M R2 ,M Ri ), (17) 

with 

M Rl = (A 2 v 16 ) 2 /M 2 , M R2 = (A 2 v^f/Mx, M R:i = (A 3 v 16 ) 2 /M 3 . (18) 



3 Procedure 

Renormalization Group Equations 

The model parameters, summarized in Tab. 1, are defined at the grand unification scale Mq with 
the exception of tan (3 and pi that are defined at the electroweak scale. At the GUT scale, «g = 
ai(Mc) = a 2 (Mc) and as(Mc) = «g(1 + £3), where £3 is the GUT scale threshold correction 3 
necessary to fit the strong coupling to experimental data at the electroweak scale, Mz- These 3 
gauge parameters, the 11 Yukawa textures (described in Section 2.2), 5 SUSY boundary conditions, 
and 3 real neutrino mass parameters allow us to completely define the model at the GUT scale and 
derive all the parameters of the minimal supersymmetric standard model (MSSM). 

3 Without presenting a complete GUT we leave £3 as a free parameter. In this way, our analysis will also apply to 
orbifold GUTs or string compactifications with a scale of order M G - 
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First, the GUT scale parameters are RGE evolved to the right-handed neutrino scale where the 
RH neutrinos are integrated out (see Fig. 2 on the following page). The right-handed neutrinos 
have three different scales associated with them, and the most relevant one is the third-family RHN 
that is mostly responsible for splitting the up and down type Higgs masses. We therefore choose to 
integrate out all the right-handed neutrinos at one single scale, M/v r = Mjj 3 . 

Below the scale of the RHNs, we use the 2-loop MSSM RGEs for both dimensionful and dimension- 
less parameters. Ideally, one should evolve all parameters to the scale of the heavy scalars (mie in 
this case, as shown in Fig. 2) and integrate them out and proceed to evolve to the weak scale using 
an effective theory without the first two generation scalars. We choose an alternative approach 
and use the 2-loop MSSM RGE 4 evolution down to the weak scale and correct for the additional 
running by including 1-loop threshold corrections to the relevant observables 5 . This approximation 
eliminates the need to define multiple effective theories. In our analysis, we have been careful to 
take into account the corresponding threshold corrections for all observables. 

Electroweak Observables 

At the weak scale, we calculate the SUSY spectrum and the SUSY threshold corrections to the 
fermion masses and CKM matrix elements. Especially in the large tan /3 regime, these SUSY 
threshold corrections are very important for the down type quarks and charged leptons and can 
be at the percent level in Yukawa-unified SUSY models [21]. We then use the threshold corrected 
fermion masses to determine the tree level masses for the squarks and sleptons. In addition, we also 
determine the one-loop pole mass for gluino and the CP-odd Higgs mass. The precision electroweak 
observables Mz, Mw, G^, a~^(Mz), a s (Mz) are calculated including 1-loop threshold corrections, 
using the procedure described in [35, 36]. Following the prescription in [35], the condition for 
consistent radiative electroweak symmetry breaking is also imposed at the weak scale, and for 
this, we use the physical Z pole mass. The parameter \i is fixed by this procedure via a separate 
X 2 minimization, and in the process, we fit the Z mass precisely to the physical Z pole mass. 
In the calculation of Mz and Mw, we only include the 1-loop corrections from the third family 
scalars, since the first two generation scalars are integrated out at mie- We assign a theoretical 
uncertainty of 0.5% to our calculation of the electroweak observables (except for Mz) due to the 
approximate treatment of thresholds described above. We also assign a 1% theoretical uncertainty 
to our calculation of G M , since we neglect the SUSY vertex and box diagrams. Finally, to compare 
to experiment, a em is evolved to zero momentum transfer. 

Charged Fermion masses and mixing angles 

Below Mz , we integrate out all SUSY partners and electroweak gauge bosons to obtain an effective 
SU(3) x U(l) om low energy theory. We use 1-loop QED and 3-loop QCD RGEs to renormalize to 
the appropriate scales and calculate the low energy observables. We fit the top quark pole mass, 

4 In scenarios with heavy scalars, it has been shown that the 2-loop contributions to the third generation scalars 
can lead to dramatic consequences, like driving the stop mass squared negative [27] and thus it is important to include 
the 2-loop RGEs in scenarios such as discussed here. 

5 For the calculation of Higgs mass, we define an effective theory at the scale Msusy and interface our calculation 
with the code by authors in Ref.[28] 
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_GUT (24 parameters) 

Model defined in terms of 24 real parameters: ac, Mg, £3, wii6, 
M1/2, A , ra Hu , m Hd , A, e, e' , p, a, e, f, M Rl , Mr 2 , Mr 3 , tan/3, /i 



RGEs for MSSM w/right-handed neutrinos 



-RHN 

Right-handed neutrinos integrated out 



RGEs for MSSM 



-m 16 _ 

1st and 2nd generation scalars are integrated out 



SUSY 



SUSY spectrum & flavor observ- 
ables calculated (susy_f lavor) 



/ 



RGEs for MSSM w/o 1st 
and 2nd generation scalars 



RGEs for MSSM w/o 1st 
and 2nd generation scalars 



Corrections 




^M top 

Calculate top pole mass 



RGEs for SM 



Mi 



EW 



SM spectrum & flavor observables cal- 
culated (Superlso) & tan/3, /i, m h , 
m H , m A , m H ± 



Compare 



RGEs: 3-loop QCD 
1-loop EW 



_2 GeV 

Calculate light quark 

masses m u , m d , m s 



• Compare ■ 



r 



Compare 

Experiment (36 Observables) . 

Gauge & EW sector: M z , M w , a em , G^, a 3 , M h 

Quark sector: M t , m b , m c , m B , m d /m s , 1/Q 2 , \V US \, \V cb \, 
IK&I, \V td \, |Vy,sin2/3 

Lepton sector: M T , M^, M e , 612, #23, #13, m^i, m 31 

Flavor observables: e K , AM Ss /AM M , AM B d, B X,7, 
B s /t + /i _ , B d -¥ B — > tv, B —¥ K*fi + fi~ (3x) 



Figure 2: This schematic shows the steps that must be employed to evolve a GUT model to the low energies and 
calculate observables at the relevant scales to compare with experimental data. Note that we use threshold corrections 
instead of integrating out two family of scalars from the particle spectrum. 



Observable 


Exp. Value 


Ref. 


Program 


Th. Error 


a 3 {M z ) 


0.1184 ±0.0007 


[29] 


mat on 


0.5% 




1/137.035999074(44) 


[291 


mat on 


0.5% 


I* 


1.16637876(7) x 10~ 5 GeV -2 


[291 

L J 


mat on 


1% 


M w 


80.385 ± 0.015 GeV 


[29] 


mat on 


0.5% 


M z 


91.1876 ±0.0021 


[29] 


Input 


0.0% 


M t 


173.5 ± 1.0 GeV 


[29] 

L J 


mat on 


0.5% 


Tflh(Tflh) 


4.18 ± 0.03 GeV 


[291 

L J 


mat on 


0.5% 




1.275 ± 0.025 GeV 


[291 


mat on 


0.5% 


m s (2 GeV) 


95 ± 5 MeV 


[291 

L J 


mat on 


0.5% 


m s /m d (2 GeV) 


17- 22 


[291 


mat on 


0.5% 


Q 


21 - 25 


[291 


mat on 


5% 


| * us\ 


0.2252 ± 0.0009 


[291 


mat on 


0.5% 


\v ub \ 


0.00377 ±0.00085 


[29] 


mat on 


0.5% 


\v cb \ 


0.04065 ±0.00195 


[29] 


mat on 


0.5% 


\v td \ 


0.00840 ± 0.0006 


[29] 


mat on 


0.5% 


\Vts\ 


0.0429 ± 0.0026 


[29] 


mat on 


0.5% 


sin 2/3 


0.679 ±0.020 


[29] 


mat on 


0.5% 


M r 


1776.82 ± 0.16 MeV 


[29] 


mat on 


0.5% 




105.6583715(35) MeV 


[29] 


mat on 


0.5% 


M e 


0.510998928(11) MeV 


[29] 


mat on 


0.5% 


M h 


125.3 ±0.4 ±0.5 GeV 


[30] 


Ref. [28] 


3 GeV 


sin 2 On 


0.27 - 0.34 (3(7 range) 


[31] 


mat on 


0.5% 


sin 2 2 3 


0.34 - 0.67 (3(7 range) 


[31] 


mat on 


0.5% 


sin 2 9i3 


0.016 - 0.030 (3(7 range) 


[31] 


mat on 


0.5% 


Am 2 i 


(7.00 - 8.09) x 10" 5 eV 2 (3(7 range) 


[31] 


mat on 


0.5% 




(2.27 - 2.69) x 10" 3 eV 2 (3cr range) 


[31] 


mat on 


0.5% 


BR(6 -> s-y) 


(343 ± 21 ± 7) x 10~ 6 


[32] 


Super I so 


(181 - 505) x 10~ 6 


BR(B — > K* (iii) i^„2<-<; r.y2 


(1.97 ±0.21) x 10~ 7 


[32] 


Super I so 


(0.79 - 3.15) x 10~ 7 


BR( _L> — } JV WW n io/ Vic u2 

V ^rV14.18<g-<lo GeV 


1.20 + n !n x 10~ 7 

— U.lU 


[321 

L J 


Super I so 


(0.48 - 1.92) x 10~ 7 


? 2 (A FB (B -> if/u/u)) 


4.9±l;| GeV 2 


[33] 


Super I so 


4.86 - 4.94 


BR(S S -> v+vT) 


3.2 x 10~ 9 


[34] 


susy_f lavor 


1.5 x 10~ 9 


BR(S„ -> ri/) 


(166 ± 33) x 1Q- 6 


[32] 


susy_f lavor 


(83 - 249) x 10~ 6 


BR(B d -> 


< 8.1 x 10~ 10 


[32] 


susy_f lavor 


< 9.72 x 10~ 10 


Am B(i 


(3.337 ±0.033) x IO" 10 MeV 


[29] 


susy_f lavor 


(2.67 - 4.00) x 10" 10 


Am Bs / Am Bi 


35.06 ±0.42 


[29] 


susy_f lavor 


28.05 - 42.07 




(2.228 ±0.11) x 10" 3 


[29] 


susy_f lavor 


(2.00 - 2.45) x 10 -3 



Table 2: The 36 observables that we fit and their experimental values. In the 4th column, we indicate the software 
package that gives us the theoretical prediction. In the last column, we show what we have assumed for the theoretical 
errors. Here, Q 2 — (m 2 — l/4(m„+m,j) 2 )/(mj - m 2 ) is defined on p. 657 of Ref. [29]. The number(s) in brackets after 
some of the values indicate the la uncertainty in the last digit(s). Capital letters denote pole masses. We take LHCb 
results into account, but use the average by Ref. [32]. All experimental errors are la unless otherwise indicated. To 
account for the inconsistencies in the inclusive and exclusive measurements of \V u b\ and IVc&l, we allow our result to 
be within the experimental error from both the inclusive and the exclusive measurement. To minimize theoretical 
uncertainties, we fit the ratio Ams s /Am_B d and derive its error by the usual formula for error propagation using the 
value Am_B s = (117.0 ± 0.8) x 10 -10 MeV [29] and assuming no correlations between the errors. Finally, the Z mass 
is fit precisely via a separate \ 2 function solely imposing electroweak symmetry breaking. 
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and the bottom and charm quark MS masses are calculated at their respective masses. All the 
other light quark masses are calculated at the scale of 2 GeV. We fit 7 observables relevant to quark 
masses, 3 charged lepton masses, and 6 CKM observables. The theoretical uncertainty in their 
calculation is again estimated to be 0.5 %. Since the light quark masses are not measured to very 
high precision, we choose to fit multiple correlated observables. These include the MS strange quark 
mass, the mass ratio md/m s and the mass ratio Q defined in the PDG [29] as 



n.2 m 2 s -l A(m u + m d ) 2 . I m u \ 1 ( m s \ 

Q = o 2 ' or equivalently, +7^ — = 1 ( 19 

m z d -mi \m d J Q l \m d J 

The CKM matrix is calculated from the left and right mixing matrices by diagonalizing the Yukawa 
matrices and including the SUSY threshold corrections. 6 CKM observables (|V^ S |, \V u b\, \V c b\, 
\Vtd\j \Vts\ an d sin 2/3) are included in our global fit analysis. To account for the inconsistencies in 
the inclusive and exclusive measurements of \V u b\ and |V^,|, we allow our result to be within the 
experimental error from both the inclusive and the exclusive measurement. The pole masses in the 
lepton sector are calculated with 1-loop electromagnetic threshold corrections. 

To execute the steps elaborated so far, we use a code maton, originally developed by Radovan 
Derrmsek to study Yukawa unification in the SO(10) model with D% x [U(l) x x Z3] family 
symmetry [17]. maton has been restructured and extended appropriately to adapt to the current 
analysis. 



Higgs Mass 

The recent observation of the Higgs boson at the LHC [30, 37] will allow us to highly constrain 
the parameter space of the model. Flavor constraints have already pushed the first two generation 
scalars of Yukawa- unified SUSY models > 10 TeV [38]. In contrast, the third family scalars have 
mass about a few TeV, purely by the effects of RGE running. The hierarchy between the first two 
and the third generations alleviates the constraints from flavor physics and CP violating observables, 
and at the same time eases the large fine-tuning in models with heavy scalars. In addition to the 
TeV range scalars, the large ^4-terms make it easy to obtain a Higgs mass of about 125 GeV. We 
integrate out all the scalars (including the third generation squarks and sleptons) below the scale 
Msusy, and calculate the Higgs boson mass using the dedicated code by the authors of Ref. [28], 
that is best suited to our case where the sfermions are very heavy. Given the boundary conditions, 

fi{M z ), Mi(Msusy), M 2 (M SU sy), M 3 (M SU sy), M SU sy, tan/3, A t (Msusy) 

at the scale Msusy = a/M^ x M^ 2 , (where Mi, M2, M3 are the gaugino masses at the scale 
Msusy) > the routine [28] determines the Higgs mass by calculating the corrections to the Higgs 
quartic coupling: 

M h = ^ A( ^ SY) (1 + S SM (Msusy) + <5 x (M SUS y)) (20) 

5 X are the contributions from chargino and neutralino diagrams. The quartic coupling A (Msusy) 
is given by: 



A(M S usy) = 7 (s 2 +S /2 ) + 



9 2 +9' 2 \ Xl Xf 



8/i| / M| USY 12M| USY 



(21) 
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We have to point out an important difference in our approach. The conventional method is to use 
the SM inputs of Mz, G^, aJ^(Mz), a s (Mz), M t , m&(rri6), M T to determine the gauge and the 
Yukawa couplings at the scales Mgusy and further constrain the GUT scale parameters. We instead 
like to predict these low energy observables and constrain the GUT scale parameter space based on 
a global x 2 fit to the data. In our calculation of the Higgs mass, we take the gauge and Yukawa 
couplings as input at the scale MgusY> obtained from RGE evolution using maton and calculate the 
Higgs mass using these inputs. The approach we adopt here is purely top-down. We have adapted 
the routine [28] to suit this line of analysis. Nevertheless, we have compared the spectrum we obtain 
from maton with that from sof tsusy 6 [39] and find good agreement. 

Neutrino Sector 

We are fitting 5 observables in the neutrino sector: the mixing angles 612, #23) #13, and the mass- 
squared differences A77131 = m\ — m\ and Arn-21 = m?, — mf (cf. Tab. 2). The most dramatic change 
in the experimental determination of the neutrino parameters in recent years comes from the Daya 
Bay and Reno collaborations [40, 41] that have confirmed that #13 ~ 9° is indeed large. Moreover, 
there are tentative hints that 823 is not maximal [42, 43]. Whereas Ref. [42] sees a preference at 
~ 2<7 — 3a for the first octant, i.e. 623 < 45°, Ref. [43] finds an equal probability for 623 being larger 
or smaller than 45°. In the following, we will be using the best-fit values and the 3a uncertainties 
quoted by the NuFIT collaboration [43] which are in agreement with Ref. [42] at 3a. 

Flavor Physics 

The strongest constraints on the model come from B-physics. For calculating the flavor observables, 
we use two publicly available codes, namely susy_flavor [44] and Superlso [45, 46]. Since the 
boundary conditions that we impose at the GUT scale may generate large off-diagonal and in 
general complex entries at the low scale, susy_f lavor is better adapted to our needs. Note that 
susy_f lavor, in contrast to comparable programs that calculate similar processes, does not assume 
minimal flavor violation (MFV), and allows for general, full three family, complex soft parameters. 
This is particularly important in our case, since we are calculating several CP violating observables 
and need to take into account 7 the complex phases in the soft parameters. Hence, susy_f lavor is 
our default choice for all flavor observables with the following exceptions. For B — > X s ~f, we use 
Superlso, since susy_f lavor does not include the NNLO SM corrections. We have verified that 
the discrepancy between susy_f lavor and Superlso in the parameter space that is of interest to 
us is at most 10% and typically less than 7%. Also, we use Superlso for the observables connected 
to the decay process B — > K*fi + fi~, since susy_f lavor does not provide them. It is important 
to note that Superlso has some built-in assumptions that prove to be too restrictive in our case. 
E.g. Superlso assumes all soft parameters to be real, and only takes the diagonal entries of the 
third-family trilinear couplings into account. As a consequence, we have assigned larger theoretical 

6 Without making significant changes to softsusy or other publicly available codes, we find that we can only make 
rough comparisons of the spectra. This is because to the best of our knowledge, most of the currently available codes 
do not handle complex parameters. In addition, many do not include right-handed neutrinos, and do not offer an 
easy way to implement the particular GUT scale Yukawa texture of the model. 

7 We calculate the particle spectrum using maton, see comments on p. 9. To the best of our knowledge, there is 
currently no publicly available spectrum generator that fully takes into account all the complex phases of the MSSM. 
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uncertainties to the values calculated by Superlso (see Tab. 2). Additional sources of uncertainties 
in the flavor observables derive from the theoretical determination of the B meson decay constant 
and from the experimental measurements of the CKM matrix elements. 

LHCb has recently measured [34] the J5r(B s — > which is in good agreement with the SM 

prediction. This pushes the CP-odd Higgs mass to a few TeV and hence leads to the Higgs decoupling 
limit. Thus the light Higgs is predicted to be SM-like. The recent observation of zero-crossing in the 
forward-backward asymmetry of B — > K*fi + /Jj~ constrains the Wilson coefficient C-j to be of the 
same sign as that in the SM. This imposes the additional constraint for the model that if \x > 0, 
in order to satisfy the branching fraction observed in the process B — > X s ^f the first two generation 
scalars have to be heavier than at least 10 TeV. 

Global Fit 

In the last step of our calculation, we construct a \ 2 function in terms of the 36 calculated observables 
(see Tab. 2). 

X 2 = £ | ' (22) 

i 1 

Hi and yf ata are the theoretical prediction and experimental measurement, respectively, for each 
observable. o~i is the error on each observable, the theoretical and experimental errors added in 
quadrature. In the general case, we vary 23 parameters (see Tab. 1 and note that m\§ is fixed in all 
the analyses) in order to fit 36 observables, which amounts to 12 (or 13 counting the separate fit to 
the Z pole mass) degrees of freedom (d.o.f.). We will consider the \ 2 P er d.o.f. for the model as a 
qualitative measure of the goodness of fit. We will look at the pulls from the individual observables 
to assess the goodness of fit of the model. 

Finding the global minimum for a model with 23 parameters is a formidable task. In the present 
analysis, we minimize the \ 2 function using the Minuit package maintained by CERN [47]. Note 
that Minuit is not guaranteed to find the global minimum, but will in most cases converge on a 
local one. For that reason, we iterate 0(100) times the minimization procedure for each set of 
input parameters, and in each step we take a different initial guess for the minimum (required by 
Minuit) so that we have a fair chance of finding the true minimum. This, of course, requires large 
computing resources, and to that end we have used the Ohio Supercomputer Center in Columbus 
and the "Centre de Calcul de l'lnstitut National de Physique Nucleaire et Physique des Particules" 
in Lyon. 

4 Third family analysis 

In this section we analyze the consequences of Yukawa unification for the the third family in the 
context of minimal SO (10) supersymmetric grand unification defined by the superpotential term, 
W D A I63 10 I63. The aim of this analysis is to study the SUSY spectrum, and we argue that 
the constraints on the SUSY spectrum come predominantly from the third family, the lightest 
Higgs mass and the branching ratio BR(B S — > /i + )• There are 24 parameters in total in the 
DR model [16, 17], and in this section we focus on 11 parameters (summarized in Tab. 1) that 
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are used to evaluate 11 low energy observables, M\y, Mz, Gf, ®em> a s(Mz), M t , m&(m&), M T , 
Br(B — > Xsj), Br(B s — > fi + fi~), and the lightest Higgs mass, Mh- We specify the model with 
the full 24 parameters, but we only vary 11 in the minimization procedure to fit the 11 observables 
listed above. The irrelevant parameters for this analysis, namely, the neutrino parameters and the 
off-diagonal Yukawa textures, are set to constant values and do not enter into the minimization 
procedure. 8 Similarly, the low energy observables connected to the first two families do not enter 
the x 2 function. The effects of the off-diagonal Yukawa textures will be discussed in Section 5. 



— 3rd family analysis w/_U, > ssn GeV 

— Third family analysis wJMg > 1000 GeV 



— 3rd family analysis w/Afj > 850 GeV 

- - 3rd family analysis w/JWj > 1000 GeV 
- . - 3rd family analysis w/.W$ > 2000 GeV 
3rd family analysis w/JWj > 3000 GeV 



CM 1 1 1 1- 0-1 , , 

innnn 15000 20000 25000 30000 10000 15000 20000 25000 30000 

mi fl [GeV] mie [GeV] 

(a) With increasing mw, X 2 first dramatically de- (b) As we increase the lower bound on the gluino mass, 

creases, and after reaching a minimum around mis — we find that x 2 dramatically increases for constant mis- 
20 TeV, starts increasing again. 



Figure 3: % 2 vs - m i6 for the third family analysis. The red and blue filled circles correspond to minima of x 2 f° r 
the values of mi6 indicated on the x-axis, where we have interpolated between them using cubic splines for ease of 
inspection. The curves correspond to different lower bounds on the gluino mass. 

Consider first the SUSY spectrum in our analysis. The first and second family squarks and sleptons 
have mass of order mie, while stops, sbottoms and staus are all significantly lighter. This is the 
inverted scalar mass hierarchy which is a direct result of RG running. Nevertheless, gluinos are 
always lighter than the third family squarks and sleptons, and the lightest charginos and neutralinos 
are even lighter. Note the states X s1 and x\ are approximately degenerate. A detailed spectrum 
is given in Tab. 3. Recent results from CMS and ATLAS give lower bounds on the gluino mass. 
These bounds are given in terms of the CMSSM or simplified models. The simplified models which 
are most relevant for our analysis are those in which (a) the third family of squarks and sleptons 
are lighter than the first two, and (b) the gluino is lighter than the stops and sbottoms. In this 
case, the lower bound on the gluino mass is now of order 1 - 1.2 TeV, assuming the branching ratio 
BR(g -»• tt Xi) = 100% or BR(g -> bb Xi) = 100% [48, 49]. Although neither simplified model is 
appropriate for our model, we nevertheless impose a lower bound on the gluino pole mass in order 
to be roughly consistent with the latest LHC results. 

In Fig. 3(a) we present the best x 2 fits as a function of mie f° r two values of the lower bound that 
we impose on the gluino pole mass, i.e. 850 and 1000 GeV. We note that x 2 is relatively insensitive 
to these lower bounds on the gluino mass, although lower values of Mg are slightly favored. The 

8 In this section, when calculating flavor violating observables, we use susy_f lavor with the experimental input 
values for the light fermion masses and mixing angles. 
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10000 15000 20000 25000 30000 



mis [GeV] 

Figure 4: The contribution of Mh to the global y 2 as a function of the lower bound on the gluino mass (vertical axis) 
and the value of mie (horizontal axis). The Higgs mass is mainly responsible for the steep increase of \ 2 observed in 
Fig. 3(a). 




,11 -\ ■— , ■ , ■ , ■ h ' 850 41 } 1 ^ 1_1 [_ 

10000 15000 20000 25000 30000 lOOOO 15000 20000 25000 30000 



mia [GeV] mm [GeV] 

Figure 5: The best \ 2 fits for the third family analysis (left) and contours of constant gluino mass ( roughly 
horizontal lines) and stop masses (vertical lines) (right). Note, for larger values of mi6 > 22 TeV, the best fit gluino 
pole mass is always much larger than the lower bound imposed. 

minimum x 2 is found for mig = 20 TeV, and x 2 increases as miQ either decreases or increases. 
Features of the model like the large A-terms and large tan/3 are favorable to obtain a Higgs mass 
in the range of 122 - 127 GeV as observed at the LHC. However, the largest contribution to x 2 f° r 
lower values of m\Q comes from the Higgs mass constraint (see Fig. 4). 
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As the lower bound on the gluino mass is increased to 2 or 3 TeV, we find that x 2 dramatically 
increases (see Fig. 3(b)). Note, this is predominantly due to the constraint from the light Higgs mass 
(Fig. 4). The simple explanation for this fact is that as the gluino mass increases the magnitude 
of At at Msusy also increases, due to the infra-red fixed point. This has the effect of decreasing 
the light Higgs mass because now Xt > v6-Mstjsy which goes beyond maximal mixing. As a 
consequence, there appears to be an upper bound on the gluino mass of order 2 TeV, which makes 
gluinos inevitably observable at the LHC 14 TeV. However, as discussed earlier, the usual simplified 
models do not apply since gluinos decay with branching ratios g — >• ti x®i 2) > ^ 2) > ^ X7i 2) > 
Xn 2)' 9 X(i 234) which are all significant. 

In Fig. 5(a) we give the best x 2 fits for the third family analysis as a function of the lower bound on 
the gluino mass and the value of tuiq. In Fig. 5(b) we give the contours of constant gluino masses 
(roughly horizontal lines) and stop masses (vertical lines) 9 . Note, for values of miQ > 20 TeV, the 
best fit gluino pole mass is always much larger than the lower bound imposed. 



5 Full Three family analysis 

In this section we present the global x 2 analysis for three families including all 24 arbitrary pa- 
rameters. The x 2 function includes 36 observables. We present our results for fixed values of m\Q 
in Fig. 6 on the next page and in Tab. 5 to Tab. 9. 

In Fig. 6 we give the best x 2 fits for two different values of the lower bound on the gluino pole 
mass imposed in the analysis. Fig. 3 and Fig. 6 have similar behavior. The value of x 2 increases 
dramatically for values of mi6 < 15 TeV. For larger values of m\§ > 25 TeV the increase is much 
slower. In the three family analysis, the minimum x 2 occurs around m\% ~ 20 TeV, just as in 
the third family analysis. Moreover, the input parameters which minimize x 2 111 the third family 
analysis also minimize x 2 f° r the full three family analysis. 

In Tab. 5 to Tab. 9 we present the best x 2 fits for values of mi6 = 10, 15, 20, 25, and 30 
TeV, respectively. The best fit overall comes for mie = 20 TeV with x 2 /d.o.f ~ 2 (see Tab. 7). 
Let us just comment on a few of the initial values of the parameters for this point. We find 
«g 1/26, Mq ~ 3 x 10 16 GeV, and 63 = —1.45%. The magnitude of the Yukawa couplings are 
hierarchical. As expected, we have Aq ~ — 2miQ, fj,, M1/2 <C ?ni6 and tan/3 ~ 50. The average 

Higgs mass parameter and the relative splitting are given by mio = y 1 \ m 'jj u + m H d )/^ ~ ^6 TeV and 
Aj^ H = (m 2 H — m 2 Hu ) / (2m\ ) « 0.07, respectively. The gravitino mass for this model is expected 
to be of order the largest scalar mass, i.e. M 3 / 2 ~ mio ~ 26 TeV. We used mi cr Omegas [51] to 
calculate the relic abundance for the benchmark points considered and found J7 th = 22.2 (10 TeV), 
ft th = 0.776 (15 TeV), fi th = 70.0 (20 TeV), Cl th = 90.2 (25 TeV), th = 123 (30 TeV). This is a 
consequence of a purely bino-like LSP. In this case, a non-thermal process would be necessary to 
accommodate the observed dark matter abundance. Assuming the correct dark matter abundance, 
a bino LSP would not have been observed yet by direct detection methods, but should be observable 
by future detectors [52]. 

9 In a recent analysis [50] , the authors found an upper bound on the stop mass for good Yukawa unification. Their 
result is a consequence of the constraint /1 < 1000 GeV, in order to satisfy dark matter bounds. We do not make any 
such assumption and do not find an upper bound on the stop mass for mi6 < 30 TeV. 
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Figure 6: x 2 vs - m i6 for the full three family analysis, x 2 is very large for low values of mm < 15 TeV. The blue 
(solid) and red (dashed) lines correspond to a lower bound on the gluino mass of 850 GeV and 1000 GeV, respectively. 



Let us now focus on the fit. Consider the observables with the largest pulls. Roughly half the 
contribution to x 2 a t this point comes from just two observables, namely md/m s and sin 2/3. Our 
value of rrid/m s is larger than the experimental value, and this implies that our value of m u /md ~ 0.9 
(see Eq. (19)). We have allowed \V U b\ to range over values consistent with both exclusive and 
inclusive measurements. We find that our fit is more consistent with exclusive measurements. 
Moreover, our fit value of sin 2/3 is at the 3a lower bound allowed by the experiments. Otherwise we 
are able to fit an amazing array of experimental observables. The light Higgs mass is fit to within 
the ± 3 GeV theoretical uncertainty we have assigned. As for the neutrino mixing angle #13 we 
obtain a value closer to 6°, rather than the present experimental value of approximately 9°. This 
may be a problem, however, it has been noticed recently that flavor violating corrections to the 
Kahler potential can have a significant effect on #13 without affecting the other larger mixing angles 
[53]. Our neutrino spectrum corresponds to the normal hierarchy. Note that the two large mixing 
angles are a consequence of a hierarchy in the right-handed neutrino masses. 

In Tab. 3 we summarize the predictions for the SUSY spectrum given values of m\Q = 10, 15, 20, 25, 30 
TeV, respectively. We give the spectrum of the lightest squark, slepton and gaugino masses, and 
the CP odd Higgs mass Ma- The first and second generation squarks and sleptons all have mass 
of order m\§. Note, in order to fit the branching ratio BR{B S — > /i + with large tan/3, we have 
Ma 3> Mz- Thus we are in the decoupling limit where the light Higgs is predicted to couple to 
matter just like the Standard Model Higgs. Therefore, any deviation from this prediction would 
rule out our model. Finally, in Tab. 4 we present results for yet to be observed quantities such as 
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mi6 



X 



M A 

mr. 



m 



m^o 

X\ 

m~+ 

Xi 
Mn 



10 TeV 



49.65 



2333 

1681 

2046 

3851 

133 

260 

853 



15 TeV 20 TeV 25 TeV 



31.02 



26.58 



27.93 



3662 

2529 

2972 

5576 

134 

263 

850 



1651 

3975 

5194 

7994 

137 

279 

851 



2029 

4892 

6353 

9769 

149 

309 

910 



30 TeV 



29.48 



2036 

5914 

7660 

11620 

167 

351 

1004 



Table 3: SUSY Spectrum corresponding to the benchmark points presented in Appendix A. The first two generation 
scalars have mass of the order of mie- 





Current Limit 


10 TeV 


15 TeV 


20 TeV 


25 TeV 


30 TeV 


e EDM xlO 28 
H EDM x 10 28 
r EDM xlO 28 


< 10.5 e cm 

(-0.1 ±0.9) x 10 9 e cm 

-0.220 - 0.45 x 10 12 e cm 


-0.224 

34.6 

-2.09 


-0.0408 

6.23 

-0.394 


-0.0173 

3.04 

-0.185 


-0.0113 

1.77 

-0.109 


-0.0084 
1.20 

-0.0732 


BR(/z -> erf) x 10 12 
BR(r erf) x 10 12 
BR(r -»■ // 7 ) x 10 s 


< 2.4 e cm 

< 3.3 x 10 4 e cm 

< 4.4 e cm 


5.09 
58.8 
1.75 


1.23 
13.9 
0.498 


0.211 

2.40 

0.0837 


0.0937 

1.04 

0.0385 


0.0447 

0.502 

0.0182 


sin 5 




-0.60 


-0.87 


-0.27 


-0.42 


-0.53 



Table 4: Predictions from the full three family analysis. The dipole moments and branching ratios were calculated 
using susy_flavor. 



electric dipole moments of charged leptons, flavor violating processes such as BR(fi — > e 7) and 
the CP violating angle in the lepton sector, sin 5. The value of sin 5 is close to zero and is thus 
consistent with tentative emerging hints for 5 ~ tt [42]. We also find that the BR(ft — > e 7) may in 
fact be observable by the MEG experiment in a few years [54]. 



6 Summary and Conclusions 



We have performed a global y 2 analysis of an SO(10) SUSY GUT times a D 3 x [U(l) x Z 2 x Z 3 ] 
family symmetry. The model fits all fermion masses and mixing angles, as well as many flavor 
observables, quite well. The model has 24 arbitrary parameters which we use to fit 36 low energy 
observables. Five of these parameters include the soft SUSY breaking masses, a universal squark 
and slepton mass, m\%; a universal cubic scalar coupling, Aq; a universal gaugino mass, M 1 / 2 and 
split Higgs up and down masses, mn u , mn d - The model has gauge coupling unification and top, 
bottom, r, u T Yukawa unification at Mqut- We have analyzed the model for the third family alone 
and then for three families. We have shown that the SUSY spectrum is predominantly determined 
by fitting the third family and light Higgs masses and the branching ratio BR(B S — > fi + /x - ). In 
Tab. 5 - Tab. 9 we give the best three family fits for five different values of the universal scalar mass 
mi6. The best overall fit is found for mis ~ 20 TeV. The SUSY spectrum for these best fit points 
are given in Tab. 3. 
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Our model makes several significant predictions. 

(i) The first and second family of squarks and sleptons obtain mass of order itiiq, while the 
third family scalars are naturally much lighter. Then gluinos and the lightest chargino and 
neutralinos are always lighter than the third family squarks and sleptons. 

(ii) Due to Yukawa unification of the third family at the GUT scale we have tan f3 ~ 50. In order 
to fit the branching ratio BR{B S — > p, + we find the CP odd Higgs mass, rriA ^> Mz- Hence 
we are in the decoupling limit and the light Higgs is predicted to be Standard Model-like. 

(iii) In order to fit the light Higgs mass, we find an upper bound on the gluino mass, M g ~ 2 TeV. 
Thus gluinos should be observable at LHC14. 

(iv) No simplified model studied to date describes the relevant gluino decay branching ratios (See 
scenarios studied in [55]). Thus in order to constrain our theory we need both CMS and 
ATLAS to provide detailed bounds on the p p — > g g production cross-section times branching 
ratios for the many different two and three body decay modes, i.e. g — >■ ti X? 12 ); ^ X(i2)' 

(v) We find BR(/j, — > e 7) ~ 10 -12 — 10~ 13 for values of toiq = 15 — 25 TeV. This may soon be 
observable at MEG [54]. 

(vi) We find the CP violating parameter in the lepton sector, sin 5 ~ 0, and the neutrinos obey a 
normal hierarchy. 

(vii) Since the first two family sleptons have mass of order m\§ we are not able to fit the muon 
anomalous magnetic moment, {g — 2)^. 

(viii) Our LSP is predominantly bino and thus assuming a thermal calculation of the relic abundance, 
we find r^o too large. 

(ix) The gravitino mass is naturally of order \/2 itiiq or for mi6 = 15—25 TeV we have nriQ ~ 20—35 
TeV. Thus the the model may avoid the cosmological gravitino problem. 

There is one obvious issue with the model regarding fine-tuning. We have not performed a detailed 
analysis of fine-tuning, but a rough measure is given by A = (j^) 2 ~ 150, corresponding to a 
fine-tuning of 1/A. As this is true for most of the surviving parameter space of the MSSM, at the 
moment we do not regard this as a serious problem. The question of electroweak fine-tuning in 
Yukawa unified models was recently studied in [56] . 

Let us now consider future directions. We will evaluate the gluino decay branching ratios in our 
model in order to compare to LHC data in a future work. In addition, we want to analyze other 
boundary conditions at the GUT scale consistent with gauge and Yukawa coupling unification. 
In particular, we will consider the "DR3" scheme [12] and also non-universal gaugino masses as 
discussed in [13], and study them with the combined predictive power of family symmetries. On 
the computational front, we would like to explore other methods to tackle the problem of finding a 
global minimum in a multi-dimensional parameter space. 

10 In order for ATLAS or CMS to test our model, we would gladly provide an SLHA2 file specifying all the low 
energy parameters of any of our benchmark points. 
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A Benchmark Points 



Table 5: Initial parameters for benchmark point with m\Q = 10 TeV: 
(l/a G , M G , e 3 ) = (25.42, 2.80 x 10 16 GeV, -2.20 %), 

(A, Ae, a, Ae, p, Ae', Ae£) = (0.61, 0.031, 1.14, 0.0048, 0.071, -0.0019, 0.0038), 

($ CT , $ g , $ 5 ) = (0.517, 0.625, 4.000, 3.497) rad, 

(mie, M 1/2 , A), //(M z )) = (10000, 239, -20247, 791.13) GeV, 

{(m Hd /m 16 ) 2 , (m H Jm 16 ) 2 , tan/3) = (1.95, 1.61, 49.42) 

(Mr,, Mr 2 , M fll ) = (3.2 x 10 13 GeV, 5.6 x 10 11 GeV, 0.9 x 10 10 GeV) 



Observable 


Fit value 


Exp value 


Pull 


a 


M z 


91.1876 


91.1876 


0.0000 


0.4559 


M w 


80.5581 


80.3850 


0.4305 


0.4022 


l/aem 


136.3909 


137.0360 


0.9415 


0.6852 


G M X 10 5 


1.1754 


1.1664 


0.7722 


0.0117 


a 3 


0.1184 


0.1184 


0.0342 


0.0009 


M t 


173.9306 


173.5000 


0.3253 


1.3238 


m b (m b ) 


4.1719 


4.1800 


0.2213 


0.0366 


M T 


1.7796 


1.7768 


0.3104 


0.0089 


m c (m c ) 


1.2782 


1.2750 


0.1246 


0.0258 




0.0962 


0.0950 


0.2437 


0.0050 


md/m s 


0.0710 


0.0526 


3.3115 


0.0055 


i/Q 2 


0.0019 


0.0019 


0.2854 


0.0001 




0.1056 


0.1057 


0.2011 


0.0005 


M e X 10 4 


5.1137 


5.1100 


0.1450 


0.0255 


\v us \ 


0.2248 


0.2252 


0.2886 


0.0014 


\Vcb\ 


0.0424 


0.0406 


0.8748 


0.0020 


\V ub \ x 10 3 


3.3462 


3.7700 


0.4985 


0.8502 


|Vtd| x 10 3 


9.5185 


8.4000 


1.8597 


0.6015 


It/ I 


U.U414 


n ri/i on 


U.oooo 


U.UVZo 


sin2/3 


0.6357 


0.6790 


2.1346 


0.0203 




0.0023 


0.0022 


0.2568 


0.0002 


AM Bs /AM B d 


59.6805 


35.0600 


3.5049 


7.0246 


AM Bd x 10 13 


3.5432 


3.3370 


0.3086 


0.6682 


m^j X 10 5 


7.6408 


7.5450 


0.1754 


0.5463 


m\ x x 10 3 


2.6521 


2.4800 


0.8182 


0.2104 


sin 2 6>i2 


0.3297 


0.3050 


0.7041 


0.0350 


sin 2 6*23 


0.6441 


0.5050 


0.8427 


0.1650 


sin 2 6>i3 


0.0128 


0.0230 


1.4585 


0.0070 


M h 


116.94 


125.30 


2.7265 


3.0676 


BR(B -» X s7 ) x 10 4 


3.9408 


3.4300 


0.3120 


1.6374 


Si?(B s -» Ai + /« _ ) x 10 9 


2.3710 


3.2000 


0.5083 


1.6308 


Bil(-B d -» n + pT) x 10 10 


1.7509 


8.1000 


0.0000 


5.2559 


572(5 -» ri/) x 10 5 


7.1988 


16.6000 


1.0525 


8.9320 


55(5 -> if>+^-)(low) xlO 8 


5.4370 


19.7000 


1.1881 


12.0051 


55(5 -» K>+^-)(high) xlO 8 


7.8844 


12.0000 


0.5651 


7.2835 


q 2 (5 -> K*fj, + fj,~) 


4.8731 


4.9000 


0.0206 


1.3009 


Total x 2 






49.6463 
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Table 6: Initial parameters for benchmark point with m\Q = 15 TeV: 
(l/a G , M G , e 3 ) = (25.50, 2.96 x 10 16 GeV, -2.40 %), 

(A, Ae, a, Ae, p, Ae', Ae£) = (0.61, 0.031, 1.14, 0.0049, 0.070, -0.0019, 0.0037), 

$g, $ p , = (0.527, 0.635, 3.881, 3.429) rad, 
(mie, M 1/2 , A), M(Mz)) = (15000, 201, -30639, 513.07) GeV, 
((m^ d /m 16 ) 2 , (m H „/mi 6 ) 2 , tan/3) = (1.97, 1.62, 49.59) 
(Mrj, Mr 2 , M fll ) = (4.2 x 10 13 GeV, 6.1 x 10 11 GeV, 1.0 x 10 10 GeV) 



Observable 


Fit value 


Exp value 


Pull 


Sigma 






91 1876 


91 1876 


n nnnn 


D 4559 

\J .'-trJOcf 






80 5671 


80 3850 


4527 


0.4022 


1 /rv 




136 4172 


1 37 0360 

1U 1 . UuUU 


9031 


6852 


Ci y in 5 




1 1 766 


1 1 664 


8739 


01 1 7 






0.1185 


0.1184 


0.1342 


0009 

yj . yj \j \j u 






173 5253 

J- f ' • -7 — -7 '.J 


173 5000 

A- 1 f J • 7 V / V 7 V ' 


0.0191 


1.3238 






4 1 QD3 


4 1 800 


D 981 3 


D 0366 






1.7756 


1.7768 


0.1366 


0089 

v 7 • V 7 v 1 \ 7 t J 


m c \m c ) 




1.2613 


1.2750 


0.5312 


0258 

V 7 • V 7 7 V .) 






D 0964 


\J. VJcfOyJ 


D 9766 


n nn5n 


iii'd/ rn s 

i /n 2 




0686 


0526 

\J .\J'J £j\J 


2 8819 


0055 




n DDI 8 

\J . \J\J ± o 


D DDI 9 


n 49nn 


n nnm 


M, 




1056 


O 1057 

yj . iu<j t 


0748 


0005 


M e X 10 4 




5.1135 


5.1100 


0.1386 


0.0255 


\v us \ 




0.2243 


0.2252 


0.6542 


0.0014 


\Vcb\ 




0.0410 


0.0406 


0.1681 


0020 


\V ub \x 10 3 




3.1115 


3.7700 


0.7745 


8502 


\Vtd\ x 10 3 




8.8886 


8.4000 


0.8124 


0.6015 






0.0401 


0.0429 


1.0638 


0026 


sin2/3 




6220 


0.6790 


2.8094 


0203 






0023 


0.0022 


0.1079 


0002 

V 7 • V 7 V 7 V t 


AM Bs /AM Bd 




37.6694 


35.0600 


0.3715 


7.0246 


AM Bd x 10 13 




4.0059 


3.3370 


1.0010 


0.6682 


m\ x x 10 5 




7.5155 


7.5450 


0.0540 


0.5463 


mgj x 10 3 




2.5097 


2.4800 


0.1413 


0.2104 


sin 2 6>i2 




0.2994 


0.3050 


0.1600 


0.0350 


sin 2 #23 




0.7414 


0.5050 


1.4323 


0.1650 


sin 2 0i3 




0.0147 


0.0230 


1.1908 


0.0070 


M h 


122.21 


125.30 


1.0080 


3.0676 


BR{B -> V s7 ) 


x 10 4 


3.5456 


3.4300 


0.0706 


1.6374 


/3i?(i3 s -»• 


") x 10 9 


4.3688 


3.2000 


0.7167 


1.6308 


/3i?(i3 d -> 


-) x 10 10 


1.3486 


8.1000 


0.0000 


5.2559 


Si?(5 -»• ri/) x 10 5 


6.2875 


16.6000 


1.1546 


8.9320 




M~)(low) xlO 8 


5.0499 


19.7000 


1.2203 


12.0051 


5i2(.B -»• iY>+ 


At - ) (high) xlO 8 


7.5449 


12.0000 


0.6117 


7.2835 


-) 


4.5922 


4.9000 


0.2366 


1.3009 


Total x 2 








31.0266 
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Table 7: Initial parameters for benchmark point with mi6 = 20 TeV: 
(l/a G , M G , e 3 ) = (25.90, 3.13 x 10 16 GeV, -1.45 %), 

(A, Ae, a, Ae, p, Ae', Ae£) = (0.60, 0.031, 1.14, 0.0049, 0.070, -0.0019, 0.0038), 

$g, $ p , = (0.533, 0.548, 3.936, 3.508) rad, 
(mie, M 1/2 , A), M(Mz)) = (20000, 168, -41087, 1163.25) GeV, 
((m^/m 16 ) 2 , (m H „/mi 6 ) 2 , tan/3) = (1.85, 1.61, 49.82) 
(Mrj, Mr 2 , M fll ) = (3.2 x 10 13 GeV, 6.1 x 10 11 GeV, 0.9 x 10 10 GeV) 



Observable 


Fit value 


Exp value 


Pull 


Sigma 






91 1876 


91 1876 


0000 


4559 






80 5452 


80 3850 


3982 


0.4022 


1 /rv 




1 37 0725 


1 37 0360 


0533 


6852 


Ci y in 5 




1 1 71 3 


1 1664 


4250 


01 1 7 






0.1184 


0.1184 


0.0467 


0009 






174.0184 


173 5000 


0.3916 


1.3238 






4 1 84Q 


4 1 800 


1 334 


0366 






1.7755 


1.7768 


0.1462 


0089 


m c \m c ) 




1.2547 


1.2750 


0.7876 


0258 






0964 


0950 


O 9807 


0050 


iii'd/ rn s 

i /n 2 




0692 


0526 


2 9891 


0055 




001 8 


001 9 


4749 


0001 


M, 




1056 


1 057 


O 1 049 


0005 


M e X 10 4 




5.1122 


5.1100 


0.0862 


0.0255 


\v us \ 




0.2243 


0.2252 


0.5964 


0.0014 


\Vcb\ 




0.0415 


0.0406 


0.4511 


0020 


\V ub \x 10 3 




3.2023 


3.7700 


0.6678 


8502 


\Vtd\ x 10 3 




8.9819 


8.4000 


0.9675 


0.6015 






0.0407 


0.0429 


0.8518 


0026 

\J • \J\J jLjKJ 


sin2/3 




0.6304 


0.6790 


2 3959 


0203 






0023 


0.0022 


3823 


0002 


AM Bs /AM Bd 




39.4933 


35.0600 


0.6311 


7.0246 


AM Bd x 10 13 




3.9432 


3.3370 


0.9072 


0.6682 


m\ x x 10 5 




7.5126 


7.5450 


0.0593 


0.5463 


x 10 3 




2.4828 


2.4800 


0.0135 


0.2104 


sin 2 6»i2 




0.2949 


0.3050 


0.2880 


0.0350 


sin 2 #23 




0.5156 


0.5050 


0.0640 


0.1650 


sin 2 0i3 




0.0131 


0.0230 


1.4134 


0.0070 


M h 


124.07 


125.30 


0.4010 


3.0676 


BR{B -> V s7 ) 


x 10 4 


3.4444 


3.4300 


0.0088 


1.6374 




") x 10 9 


1.6210 


3.2000 


0.9682 


1.6308 


5i?(5 d -> 


-) x 10 10 


1.0231 


8.1000 


0.0000 


5.2559 


5i?(5 -»• ri/) x 10 5 


6.3855 


16.6000 


1.1436 


8.9320 


5^(5 -> 


M~)(low) xlO 8 


5.1468 


19.7000 


1.2123 


12.0051 


£#(£ -»• K>+ 


At - ) (high) xlO 8 


7.7469 


12.0000 


0.5839 


7.2835 


-) 


4.5168 


4.9000 


0.2945 


1.3009 


Total x 2 








26.5812 
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Table 8: Initial parameters for benchmark point with m\Q = 25 TeV: 
(l/a G , M G , e 3 ) = (25.83, 4.17 x 10 16 GeV, -2.55 %), 

(A, Ae, a, Ae, p, Ae', Ae£) = (0.61, 0.031, 1.17, 0.0049, 0.070, -0.0019, 0.0037), 

$g, $ p , = (0.513, 0.542, 3.969, 3.503) rad, 
(mie, M 1/2 , A), M(Mz)) = (25000, 158, -51365, 1348) GeV, 
((m^/m 16 ) 2 , (m H „/mi 6 ) 2 , tan/3) = (1.86, 1.61, 49.98) 
(Mrj, Mr 2 , M fll ) = (3.2 x 10 13 GeV, 6.1 x 10 11 GeV, 0.9 x 10 10 GeV) 



Observable 


Fit value 


Exp. value 


Pull 


Sigma 




91 1 876 


91 1876 

al.lOlU 


O 0000 


4559 


iv i yy 


80 6192 


80 3850 


5824 


0.4022 


1 / ry 

L / "em 


1 37 1 624 


1 37 0360 


1844 


6852 


(7 x 10 5 


1 1 754 

1.11 (J^fc 


1 1664 


7749 


01 1 7 


/To 


0.1185 


0.1184 


0889 


0009 


iv± t 


174.5241 


173 5000 


0.7735 


1.3238 


> lb b\ ul >b) 


4 1 78Q 


4 1 800 


O 0307 


0366 


M T 


1.7761 


1.7768 


0800 


0089 




1.2529 


1.2750 


8559 


0258 


llbg 


0963 


O 0950 


9659 


0050 


m 7 / 9Ti 
" L dl " t s 

1/ V 


0702 


0526 


3 1 726 

(J.J- J \) 


0055 


001 Q 


O OOI 9 


3379 


0001 




1 057 


1057 


1 533 

VJ . 1 tJtJ*J 


0005 


M e X 10 4 


5.1102 


5.1100 


0.0083 


0.0255 


ivj 

\ ' US \ 


0.2244 


0.2252 


0.5407 


0.0014 


1 K c6| 


0.0411 


0.0406 


2090 

VJ . VJ O \J 


0020 


iKhl x 10 3 

I v | 


3.1806 


3.7700 


6933 


8502 


Vwl x 10 3 

I v ta\ XVJ 


8.9530 


8.4000 


0.9193 


0.6015 


IV 1 


0.0402 


0.0429 


1.0358 


0026 


sin2/3 


0.6318 


0.6790 


2 3268 


0203 


e# 


0.0024 


0.0022 


8902 

V / • V.- 1 ' / V 7 — 


0002 




35.1576 


35.0600 


0.0139 


7.0246 


AMpj x 10 13 


4.1075 


3.3370 


1.1531 


0.6682 


m\ x x 10 5 


7.5325 


7.5450 


0.0229 


0.5463 


x 10 3 


2.4814 


2.4800 


0.0066 


0.2104 


sin 2 6»i2 


0.2978 


0.3050 


0.2069 


0.0350 


sin 2 023 


0.5109 


0.5050 


0.0358 


0.1650 


sin 2 6»i3 


0.0140 


0.0230 


1.2789 


0.0070 




125.21 


125.30 


0.0293 


3.0676 


RR(£ -> X s ~f) x 10 4 


3.4074 


3.4300 


0.0138 


1.6374 


-»• fi + fi-) x 10 9 


2.6112 


3.2000 


0.3610 


1.6308 


B12(.B d -»• n+fi-) x 10 10 


1.0779 


8.1000 


0.0000 


5.2559 


£#(£ -»• TV) X 10 5 


6.4123 


16.6000 


1.1406 


8.9320 


BR(B K*^ + Li-) (low) xlO 8 


5.0511 


19.7000 


1.2202 


12.0051 


-»• A'* ! u + ! u~)(high) xlO 8 


7.6223 


12.0000 


0.6010 


7.2835 


4.4839 


4.9000 


0.3198 


1.3009 


Total x 2 






27.9288 
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Table 9: Initial parameters for benchmark point with m\Q = 30 TeV: 
(1/q g , M g , e 3 ) = (25.86, 4.36 x 10 16 GeV, -2.81 %), 

(A, Ae, a, Ae, p, Ae', Ae£) = (0.62, 0.031, 1.18, 0.0050, 0.069, -0.0020, 0.0037), 

($ CT , $g, $ 5 ) = (0.507, 0.534, 4.005, 3.514) rad, 

(mi 6 , M 1/2 , A , /u(M z )) = (30000, 161, -61640, 1647) GeV, 

{(m Hd /m w ) 2 , (m H Jm 1(i ) 2 , tan/3) = (1.86, 1.63, 50.15) 

(AfRg, Mr 2 , M Rl ) = (3.2 x 10 13 GeV, 6.4 x 10 11 GeV, 0.9 x 10 10 GeV) 



Observable 


Fit value 


Exp. value 


Pull 


Sigma 


Mr, 


Q1 1 876 

cf ± . J- O ( VJ 


91 1876 

al.lOlU 


O 0000 

VJ .UVJVJVJ 


4559 


IV 1 y[/ 


80 6519 

(JU. VJtJ 1 U 


80 3850 


6637 

VJ . VJVJtJ i 


0.4022 


1 In 

L / "em 


1 S7 0422 


1 37 0360 

-L KJ \ • VjKJVjVj 


0091 

VJ .VJVJ cf -L 


6852 

VJ . *JO ?JZj 


(1 x 10 5 


1 1 785 


1 1664 


1 0410 


01 1 7 

VJ . VJ J. J. 1 


/To 


0.1187 


0.1184 


2850 

VJ . Zj VJ *J VJ 


0009 

vj . vjvjvj u 




175 0383 


173 5000 


1.1620 


1.3238 


> lb b\ ul >b) 


4 1 789 


4 1 800 


0484 


0S66 

VJ .VJkJXJVJ 


M T 


1.7764 


1.7768 


0.0419 


0089 




1.2504 


1.2750 


0.9540 


0258 




0979 

VJ .VJCf 1 Li 


0950 

VJ .WJUU 


4^,69 


0050 

vj . yjvjovj 


7T7 7 / TY1 

" L dl " t S 

1/ V 


071 2 

VJ . VJ % X _■ 


0526 

\J .VJkj^jVJ 


3 3536 

' ) . fj 'J (J V) 


0055 

VJ . \JVJkjKJ 


001 9 

VJ . VJVJ ± cf 


001 9 

vj . uvj ± cf 


S099 

VJ .OVJZj Cf 


0001 

VJ .VJVJVJ ± 




1 058 


1057 

vj . i vj 'j i 


1 928 

VJ . 1 iJZj o 


0005 

vj . \jvjvjtj 


M e x 10 4 


5.1097 


5.1100 


0.0125 


0.0255 


IK«I 

| ' US \ 


0.2244 


0.2252 


0.5324 


0.0014 


\ v cb\ 


0.0405 


0.0406 


0988 

VJ . tj t/ (_)(J 


0020 

vj • \jvj £i\j 


\V„h\ X 10 3 
I v uo \ 


3.1793 


3.7700 


0.6947 


8502 


\Vtj\ x 10 3 

1 v ta\ ^ xv7 


8.9071 


8.4000 


0.8431 


0.6015 


1 v ts\ 


0396 

\J . VJtJC/VJ 


0429 


1.2682 


0026 

VJ . VJVJ^iVJ 


sin2/3 


6380 


0.6790 


2.0207 


0203 


e# 


0.0022 


0.0022 


0.0877 


0002 


1 —*- ±yj - t>S 1 Dtt 


34.0021 


35.0600 


0.1506 


7.0246 


AMpj x 10 13 


4.1018 


3.3370 


1.1445 


0.6682 


m\ x x 10 5 


7.5705 


7.5450 


0.0467 


0.5463 


x 10 3 


2.4783 


2.4800 


0.0081 


0.2104 


sin 2 6»i2 


0.3057 


0.3050 


0.0213 


0.0350 


sin 2 023 


0.5036 


0.5050 


0.0087 


0.1650 


sin 2 6»i3 


0.0130 


0.0230 


1.4280 


0.0070 




125.88 


125.30 


0.1876 


3.0676 


RR(£ -> X s7 ) x 10 4 


3.3931 


3.4300 


0.0225 


1.6374 


BR(B S -»• M+/" - ) x 10 9 


2.6139 


3.2000 


0.3594 


1.6308 


B12(.B d -»• /i+^i-) x 10 10 


1.0748 


8.1000 


0.0000 


5.2559 


£#(£ -»• TI/) X 10 5 


6.4081 


16.6000 


1.1411 


8.9320 


BR(B -> K*fi + n-)(lcm) xlO 8 


4.9279 


19.7000 


1.2305 


12.0051 


£#(£ -»• K*n + n~) (high) xlO 8 


7.4423 


12.0000 


0.6257 


7.2835 


4.4707 


4.9000 


0.3300 


1.3009 


Total x 2 






29.4783 
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